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ABSTRACT 



Context. Supergranulation is a pattern of the velocity field at the surface of the Sun, which has been known about for more than fifty years, 
however, no satisfactory explanation of its origin has been proposed. 

Aims. New observational constraints are therefore needed to guide theoretical approaches which hesitate between scenarios that either invoke 
a large-scale instability of the surface turbulent convection or a direct forcing by buoyancy. 

Methods. Using the 14-Mpixel CALAS camera at the Pic-du-Midi observatory, we obtained a 7.5 h-long sequence of high resolution images 
with unprecedented field size. Tracking granules, we have determined the velocity field at the Sun's surface in great detail from a scale of 
2.5 Mm up to 250 Mm. 

Results. The kinetic energy density spectrum shows that supergranulation peaks at 36 Mm and spans on scales ranging between 20 Mm and 
75 Mm. The decrease of supergranular flows in the small scales is close to a tc 1 -power law, steeper than the equipartition Kolmogorov one. The 
probability distribution function of the divergence field shows the signature of intermittency of the supergranulation and thus its turbulent nature. 
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1. Introduction 



Supergranulation was discovered by Hart ( 1954) using Doppler 
images of the Sun. It appeared as an essentially horizontal 
flow field at a typical scale of 30 Mm. The origin of this flow 
field was first thought to be related to the second ionization 
of helium, which provides some latent heat at a depth around 
10 M m, compatible with their typical size (iSimon & Weiss 
1968J). This scenario has been much debated because of the 
weakness of the effect and the apparent vigour of the super- 
granular flow. Other ways of generating this velocity scale 
rely on large-scale instabilities of the surface turbulent flow 
( Rieutord et al . 2000), triggered by the strong density stratifica- 
tion. In this scenario, kinetic energy of granules, the small-scale 
convective cells, is thought to be p iped to a larger scale by an 
AKA-like effect (IGama et al.lll994l) . However, other kinds of 
large-scale instabilities are still possible, like a convective in- 
stability triggered by fixed-flux boundary conditions imposed 
by th e small-scale granular convection (IRincon & Rieutord 



2003). In this approach, the cooling resulting from the granu- 



lation does not suppress the convective instability of the larger 
scales; because all the heat flux is carried by the small scale, 
none is carried by the large scale instability, which shows lit- 
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tie temperature fluctuations (see lMeunier et al. 2007a). Besides 
these scenarios, recent observations of iGizon et alj d2003l) 
added some wavelike properties to supergranulation. 

Obviously, the supergranulation theory needs guidance 
from every available observational constraint. This letter 
presents the results of the observations issued from the CALAS 
project (a CAmera for the LArge scales of the Sun; see 
Meuni er et al. I l2005l) . collecting a sequence of high resolution 
wide field images, which have allowed us to capture the evolu- 
tion of a hundred supergranules at the disc centre during 7.5 h. 
We can thus give new constraints on the dynamics of the solar 
surface in the supergranulation range. After a brief description 
of the data set (Sect. 2), we show the spectral side of the super- 
granular flow (Sect. 3) as well as the probability distribution 
functions of the divergences (Sect. 4); first conclusions follow. 



2. Observational techniques and data reduction 

2.1. Data set 

On 13 March 2007, we observed the Sun at disc centre during 
7.5 h using the Lunette Jean Rosch at Pic du Midi, a 50 cm- 
refractor. Images were taken at A = 575+5nm, with a 14Mpixel 
CMOS-camera (4560x3048 pixels), with 0.115 arcsec/pixel, 
thus covering of 524x350 arcsec 2 (see Fig.[TJ. 10811 images 
were obtained with a regular cadence of one every 2.5s. Two in- 
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Fig. 1. Bottom right: the CALAS field of view on the Sun; the large rectangle indicates the size of the images, the smaller one 
shows the field where the velocities could be computed, while the small square shows the field of view of the SOT instrument on 
the Hinode satellite when tracking granules. Top: the supergranulation velocity field with the divergence contours superimposed 
(scales shorter than 8 Mm have been filtered out). A time-window of 150 min was used. Bottom left: a zoom on granulation 
showing the relative sizes of granules and supergranules. 



dependent series of ~ 1400 images were then extracted. They Because of tracking difficulties, the common field of each 
sample the same solar signal with a period of 20s but are noised series was reduced to ~400x300 arcsec 2 , thus covering a sur- 
differently by the Earth atmospheric perturbations. The com- face of 290x216 Mm 2 on the Sun (Fig.[B. 
parison between the outputs of both series allows us to evaluate 
the influence of the seeing and test the robustness of the results 
with respect to this noise. 

After recentering, subimages were k - a> filtered (with a 
threshold of 7 km/s), so as to remove, as much as possible, 
Earth atmospheric distortion , which is the main s ource of noise 
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Fig. 2. Kinetic energy spectra obtained for various time win- 
dows. The vertical dotted line indicates the position of the peak 
at 36.4 Mm. The vertical dashed line emphasizes the 10 Mm 
scale, usually taken as the upper limit of mesogranular scale. 
Two power laws are shown on each side of the peak, as well as 
the one of small-scale noise. 

2.2. Velocity fields 

Horizontal velocity fields have been obtained using the CST 
granu le tracking alg orithm dRoudier et all 19991: Rieutord et al. 
20071) . As shown in iRieutord et al.l d200ll) . granules' motions 
trace large-scale velocity fields when the scale is larger than 
2.5 Mm. Hence, we sampled the velocity field with a bin of 12 
pixels (~1 Mm). 

The velocities are obtained by tracking the granules during 
a given time window. Thus, we have access to average veloc- 
ity components, namely v x (i, j) and v y (i, j), where the overline 
refers to the time averaging imposed by the time window. This 
averaging improves the signal-to-noise ratio, but naturally de- 
creases the time resolution. Typically, the shortest time window 
that can be used is 30 min. 

Although the field is large, projection effects of the spheri- 
cal Sun are still of weak influence. At most, in the field corners, 
the correction on the velocity would be less than 4%, which is 
much less than the noise. 

We show in Fig. Q] an example of these velocity fields. The 
small-scale components of the flow (with scale s below 8 Mm) 
were filtered out using Daubechies wavelets (see Rie utord et al. 
20071) . Figure [T] shows that robust steady supergranules live 
among a wide variety of flow structures illustrating the turbu- 
lent nature of these scales and their wide spectral range. 

3. The kinetic energy spectrum 

A convenient way to view the dynamics of a flow is to ex- 
amine the spectral content of the velocity field. We therefore 



Fig. 3. Same as in Fig. [2] but for the horizontal divergence of 
the flow, d x v x + dyVy (solid line), and for the vertical vorticity, 
d x v y - d y v x (dashed line); the supergranulation peak is clearly 
visible at A — 36 Mm in the divergence spectrum but absent in 
the vorticity one. A k 2 -power law is given for comparison. 

computed the spectral density of kinetic energy E(k) associ- 
ated with the horizontal flows that we can measure. It is such 
that 
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E(k), w hich is computed in the same way as in IRieutord et al 



( 2000), is displayed in Fig.|2]for the flows determined with var- 
ious time-windows. The short time-windows of 45 min give us 
10 independent spectra, which are averaged together, while the 
whole series 450 min time-window gives us only one spectrum. 
Note that no spatial filtering was applied to the velocity field. 
We see that E(k) oc k at small scale, which is the signature of 
decorrelated random noise. 

These spectra clearly show the emergence of the spec- 
tral range of the supergranulation as the length of the time- 
window increases. Let us point out the remarkable stability of 
the wavenumber of the spectral peak when the time averaging is 
changed. This demonstrates that supergranulation is a genuine 
velocity field at the Sun's surface. It is not the consequence 
of time-averaging the small, fast turning-over small scales. Of 
course, if the averaging interval is long enough (i.e. of the or- 
der of the turn-overtime scale of supergranulation), the spectral 
peak will move to ever larger scales before disappearing. The 
maximum of the spectral density is at a wavelength of 36 Mm. 
The FWHM of the peak indicates that supergranulation occu- 
pies the range of scales of [20, 75] Mm. 

The present mean value of the diameter of supergranules is 
slightly higher than the previous determinations. For instance, 
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Meunier et al.l (12007bl) find a mean diameter of 3 1 .4 Mm with a 
technique based on the segmentation of the divergence field de- 
rived from velocities issued from a local corre lation technique 



applie d to SOHO/MDI white light images. iDel Moro et al 



(2004) also use a divergence field but derived from time- 
distance helioseismology; they find a mean size of 27 Mm, sig- 
nificantly smaller than ours. 

Both of these results are based on the divergence field and 
its segmentation. The histogram of sizes is then used to deter- 
mine the mean diameter of a supergranule. Such a technique 
necessarily underestimates the actual scale of supergranular 
flows as only part of it (the positive divergences) is used. Our 
spectra, which directly result from the measured horizontal ve- 
locities, incorporate all the components of the flow at super- 
granulation scale, and thus better reflect the dynamical state. 
Note that a blind use of the divergence spectrum in Fig.[3]would 
point to a scale of ~ 23 Mm. 

In Fig. |2] we also indicate the best-fit power laws, which 
mimic the sides of the supergranulation peak. We find that 
E(k) ~ k 3 on the large-scale side and E(k) ~ k~ 2 on the 
small-scale one. This latter power law is steeper than the -5/3 
Kolmogorov one and may be an effect of density stratification. 

To complete the spectral picture, we also show, in Fig. [3] 
the spectra of the horizontal divergence of the velocity field and 
the vertical vorticity. Beyond the A: 2 -dependence of the spectral 
densities, which is a consequence of the derivative of uncor- 
rected noise, we can see that the divergence spectrum clearly 
shows the supergranulation peak, while the vertical vorticity 
shows no signal in this range. The vorticity does show some 
weak signal however, but in the mesogranulation range, below 
10 Mm. 

4. Probability density functions 

Another way to look at the velocity fields is to consider the 
probability density functions. Such distributions have been 
computed for the velocity field, its divergence and its curl. We 
show in Fig. [4] the distribution of the divergence field when 
scales shorter than 8 Mm have been filtered out of the velocity 
field. This distribution clearly shows a positive wing larger than 
the Gaussian one, associated with a skewness of 0.38. We also 
note that the noise of the Earth atmosphere is weak enough, and 
does not perturb the result. This trend to a n exponential distri- 
bution is a signature of intermittency (e.g. Frisch|[l995l) . This 
result, with the same asymmetry between positive and nega- 
tive values, was also observed by IMeunier et al. (2007b) with 
a completely different set of data and method. Thus, the in- 
termittency of supergranulation seems to be a robust property. 
Exponential wings are usual ly less visible on the velocity (e.g. 
Vincent & Meneguzzill991 ), and, indeed, are barely noticeable 
in our data. As far as vorticity is concerned, the noise is unfor- 
tunately too high to give convincing measurements. 

5. Conclusions 

For the first time, it has been possible to follow the motion 
of well-resolved granules in a large field of view. The spec- 
tral peak of supergranulation has thus been determined in a 
very direct manner. According to our data, supergranulation 
has the most energetic motions at a scale of 36 Mm and en- 
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Fig. 4. Histograms of the divergences; solid and dashed lines 
have the same solar signal but a different noise from the Earth 
atmosphere. Scales below 8 Mm have been filtered out. A 
Gaussian distribution with the same standard deviation is over- 
plotted for comparison. 



compasses all the scales ranging from 20 to 75 Mm, as shown 
by the FWHM of the peak. Except for its amplitude, this peak is 
not sensitive to the time window used to measure the granules' 
motions. The signal also clearly appears in the divergence spec- 
tral density, but not in the vorticity; it is likely that vorticity at 
supergranulation scale near the Sun's equator is much weaker 
and does not emerge from the noise. 

Finally, our data confirm the fact that supergranulation has 
a noticeable degree of intermittency, clearly appearing in the 
distribution of positive divergence values. 

As far as the origin of supergranulation is concerned, the 
scenarios mentioned in the introduction can be tested with 
these data, either in the real space with velocity fields like the 
one shown in Fig.[T]or in the spectral space with the given spec- 
tra. 

Further work on the observational side will focus on the 
determination of the third component of the velocity field, the 
increase of the field size and the reduction of the noise, so as to 
further constrain the dynamics of scales from the granulation 
one to the 100 Mm one. 
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